High-precision metrology of highly charged ions via relativistic resonance fluorescence 
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Resonance fluorescence of laser-driven highly charged ions is studied in the relativistic regime by 
solving the time-dependent master equation in a multi-level model. Our ab initio approach based 
on the Dirac equation allows for investigating highly relativistic ions, and, consequently, provides 
a sensitive means to test correlated relativistic dynamics, bound-state quantum electrodynamic 
phenomena and nuclear effects by applying coherent light with x-ray frequencies. Atomic dipole or 
multipole moments may be determined to unprecedented accuracy by measuring the interference- 
narrowed fluorescence spectrum. 
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High-precision laser spectroscopy has resulted in cru- 
cial advancements in our understanding of nature. In 
particular, optical laser spectroscopy (OS) is a versatile 
tool to investigate correlated relativistic quantum dy- 
namics, the testing of fundamental theories like quantum 
electrodynamics (QED) (TJ [2] or parity non-conservation 
in atomic systems. The determination of atomic dipole or 
multipole moments via lifetime measurements by means 
of, e.g., visible emission spectroscopy [3], approaching the 
accuracy of one per thousand, sheds light on QED effects 
like the electron anomalous magnetic moment. Isotope 
shifts (IS) in atomic spectra which has been providing 
valuable insight into the collective structure of nuclei: 
for example, recently, isotope shifts were determined by 
collinear laser spectroscopy @]. Beyond purely nuclear 
effects, the interaction of the correlated motion of elec- 
trons and that of the nucleus can be studied in IS mea- 
surements: recently, relativistic effects on nuclear recoil 
have been measured in visible forbidden transitions of 
few-electron argon ions by a trapped-ion method 

In the regime of heavy few-electron systems, however, 
the accuracy of optical spectroscopy can seldom be ex- 
ploited due to the scarcity of low-frequency transitions. 
Therefore, one has to apply other techniques. Measuring 
x-ray emission lines of highly charged uranium ions con- 
fined in an electron beam ion trap allowed testing high- 
field QED on the two-loop level [T] and delivered a new 
value for the radius of the radioactive isotope 235 U [6]. 
Recently, a method based on the storage ring measure- 
ment of dielectronic recombination spectra yielded the 
change of charge radii for neodymium isotopes pj [5] . 

With the advent of modern short-wavelength laser sys- 
tems, the accuracy and versatility of laser spectroscopy 
may be combined with the increased sensitivity of highly 
charged ions (HCI) to relativistic and QED effects, nu- 
clear properties, and this may also give new information 
on HCI structure and dynamics relevant in astrophysi- 
cal and thermonuclear plasmas. Brilliant x-ray light has 
already enabled to study transitions in the soft x-ray 
regime involving HCI |9j. Coherent light with photon 
energies over 10 keV becomes accessible in the near fu- 



ture |10j , allowing for an extension to heavier systems and 
the exploitation of coherence properties. This would also 
ask for the validity of numerous quantum control schemes 
of resonance fluorescence [TTHT§] for concrete systems in 
the relativistic regime. 

In the present Letter, we investigate the possibility of 
measuring atomic transition dipole - or multipole - mo- 
ments and transition energies via relativistic resonance 
fluorescence of a three-level atomic configuration driven 
by two fields, namely, a short- wavelength laser and a 
long-wavelength light source in the optical regime. In 
such a three-level setting, the linewidth of the sponta- 
neous transition of interest may be rendered much nar- 
rower than the natural linewidth, with the simultaneous 
increase of the total emitted intensity by orders of mag- 
nitude. Due to this effect, the determination of atomic 
multipole moments by means of the detection of the flu- 
orescence spectrum is anticipated to increase in accuracy 
by several orders of magnitude. 

As relativistic effects on the electronic wave function 
increase rapidly with the nuclear charge number Z, we 
formulate a fully relativistic ab initio theory of coher- 
ent laser-atom interaction based on the Dirac equation. 
An approach via the time-dependent numerical solution 
of this wave equation was recently employed to describe 
ionization phenomena |20j . For our purposes, one needs 
to go beyond this approach and incorporate radiative re- 
laxation in bound-bound transitions. 

The interaction of a radiation field, having the four- 
vector potential (/x G {0,1,2,3}), with a single- 
electron atom can be described by the Hamiltonian (in 
relativistic units) Jf = J#a + + ej^A^. Here, 
and J#f are the energy operators of the atom and the 
radiation field, respectively, e is the elementary charge, 
and 7^ is the four- vector of Dirac matrices, represented in 
terms of 2 x 2 Pauli matrices. The atom is described in a 
fully relativistic manner, i.e. M'a = ij^d^ + V(r), where 
V(r) is the nuclear potential, r the distance from the 
origin, and 9 M is the space-time gradient operator. The 
Hamiltonian of the electromagnetic field J%f is given in 
terms of the photon creation and annihilation operators 
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by .fflp — u>k (^o^o-k + |^ • To put the total Hamilto- 
nian in a more convenient form, M'a and j^A^ can be ex- 
pressed in terms of the atomic operators a a b — l a )(&| de- 
scribing transitions between a state \a) and \b) (a, b e N). 
Let us introduce the coupling constants (Rabi frequen- 
cies) g a b = —fJ'ab'^'^k, proportional to the corresponding 
laser field strength, with the relativistic multipole tran- 
sition matrix element /J, ab = — {a\ct ■ £ke lk ' r |&). Here, 
ek denotes the laser polarization, u = |k| and <%k is the 
amplitude of the electric field. We evaluate the multi- 
pole transition matrix elements by using the multipole 
expansion of the vector potential 



e k e 



47T ]T 

JMX 



(Y^(fc)-e k )a^(r) 



(1) 



The electric (A = 1) and magnetic (A = 0) multipole 
potentials a^ f (r) are given by 



a j°A/( r ) = 3L(kr)YjLM(r), 
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in the transverse gauge, which can be regarded as a rel- 
ativistic generalization of the velocity form. The sym- 
bols Yj^(f) and Yjlm(?) stand for the vector spherical 
harmonics. We adopt the generalized length gauge, also 
referred to as the Babushkin gauge [21] . This general 
description also accounts for higher-order (non-dipole) 
transitions resulting in narrow spectral features. The 
Hamiltonian accounting for the interaction of the three- 
level system with two classical fields is [TTJ 
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where A x = W31— uj x and A Q = iQ 2 \ — w D are the detunings 
of the x-ray and optical laser frequencies lo x , lj q from the 
atomic transition energies. The master equation describ- 
ing the dynamics of the system with the Hamiltonian in 
the interaction picture is given by p = —i[J4fi,p] + Ap, 
where A is the Lindblad operator accounting for the re- 
laxation channels [TTJ [T7HT9] . We have also included the 
dephasing rate 7/3 of x-ray lasers by adding it to the 
decay rate 713 = |(r 31 + r 32 ) of the 3ol x-ray tran- 
sition [TTJ [22] ; since the coherence time of free electron 
lasers (FELs) is often shorter than their pulse length. r a ;, 
denotes the partial radiative rate for the a— J-b transition. 

Our goal is to evaluate the power spectrum of radiation 
scattered by a few-level atom driven by an incident field 
of arbitrary strength. To this end we introduce the field 
operators E^^r,^) for a given time t and observation 
point r in the Weisskopf-Wigner approximation 
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FIG. 1: (Color online.) Fluorescence photon spectrum for 
the 2s-f->2p3/2 transition in Li-like 209 Bi as a function of the 
fluorescence photon frequency uf. (a) Dashed (red) curve: 
An x-ray laser (I x = 5 x 10 11 W/cm 2 ) is in resonance with 
the ionic electric dipole (El) transition at U31 = 2788.1 eV 
between the hyperfine-split ground state 1 (2s with F = 4, 
Mf = 4) and the uppermost state 3 (2p 3 / 2 with F = 5, 
Mf = 5). This curve is multiplied by a factor of 5 x 10 . 
Thick (thin) dashed arrows represent fast El x-ray (slow Ml 
optical) decays, (b) Continuous (blue) curve: an additional 
optical driving (I a = 10 14 W/cm 2 ) is applied on the CJ21 = 
0.797 eV [23] Ml transition between the hyperfine-split mag- 
netic sublevels 1 (F = 4, M F = 4) and 2 (F = 5, Mf = 5). 
The inner sidebands are suppressed. See text for more details. 



and the corresponding conjugate expression for 
E(-)(r,£). Here, P(t)<+> = fi 12 (a 12 + V21) + Mi 3 (^i3 + 
CT31). The atomic dipole (or, more generally, multipole) 
is assumed to be in the plane defined by the Cartesian 
coordinates x and z, and rj is the angle between the 
dipole and the z-axis. 

We are interested in the field emitted by the atom 
fixed in position along the x-axis. The field oper- 
ator in Eq. Q can therefore be treated as scalar. 
The power spectrum S(r,uo) of the fluorescent light 
at some suitably chosen point r in the far field is ob- 
tained by taking the Fourier transform of the normally- 
ordered field correlation function (E(~) (r, t)E^ (r, t+r)) 
with respect to the time delay r. According to the 
Wiener-Khintchine theorem, the power spectrum S(u>o) 
is = iRe/ oo dr(E(-)(r,t)E(+)(r,t + T))e-^. 

We calculate the power spectrum and the transition en- 
ergies and matrix elements involved by relativistic meth- 
ods. The Dirac wave functions corresponding to a rota- 
tionally symmetric nuclear potential are given in spheri- 
cal coordinates in the form (see, e.g. [2"T] ) 
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The subscript a stands collectively for the principal quan- 
tum number n a as well as the angular momentum quan- 
tum numbers j a and l a and the magnetic quantum num- 
ber m a of a given state a = (n a ,j a , l a , ma)- For Coulomb 
potentials, the bound radial functions G a (r) and F a (r) 
are known analytically, involving confluent hypergeomet- 
ric functions with negative- integer first argument. The 
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spin-angular part of the Dirac wave function is defined 
by the spherical spinors fl Karria (f) [21 j . In the transition 
energies oj a b, also the one-loop QED corrections of self- 
energy and vacuum polarization _24] have been included. 

In Fig. [I] (a) we plot the power spectrum of resonance 
fluorescence for the case of the 2s^2p 3 / 2 electric dipole 
transition in Li-like 209 Bi (Z=83) ions. The dynamic 
(AC) Stark shift leads to a splitting of the central peak, 
giving rise to a Mollow spectrum. Due to the long lifetime 
of the upper level of the hyperfine-split ground state, level 
2, almost 100 % of the population is trapped in this level 
if only the 3f->l transition is driven coherently with an x- 
ray laser. The calculation yields for the population of the 
uppermost state a value of approx. r 2 i/(r 3 2 + 2r 2 i) ~ 
10~ 12 <C 1, resulting in a negligibly small total x-ray 
fluorescence [35] . This undesirable effect may be reversed 
if additionally the 2ol optical transition is coherently 
driven (see Fig. [IJ , leading to an efficient re-population 
of level 3. 

Furthermore, the spectral lines become substantially 
narrower due to coherence and interference effects 
(see for the pioneering non-specific treatment). The 
width of the central peak and the outer sidebands are 
given by T c = (T 31 + T 32 + j D )R + T 21 (l - R) and 
Tsb = |f(r 31 - \ 1D )R+ \T 32 {R + R 2 ) + §r 21 (i- R)\, 
respectively, with the ratio R being 331/(531 + g 2 \) [2"5] . 
This effect is shown in Fig. [I] (b) . Further increasing the 
intensity of the long- wavelength driving field and thus g 2 \ 
could even assign the narrow linewidth of 7.7- 10 -15 eV of 
the Ml hyperfine transition to the El x-ray transition of 
interest. The above line width formulas also imply that 
the dephasing width jn - typically on the order of 1 eV 
for XFELs [10] - does not hamper the observation of sub- 
natural linewidths in the x-ray regime as its contribution 
is decreased with the same factor R. 

Transition lifetimes - and related quantities like the 
atomic dipole or multipole moments - are of great in- 
terest for astrophysical applications and for testing fun- 
damental theories. Measurements of these quantities 
are particularly necessary since they are especially sen- 
sitive to the long-range behavior of atomic wave func- 
tions. Currently, even the best measurements do not 
exceed the 10 -3 level of accuracy [3J. In our scheme, the 
narrowed central and outer lines enable in principle an 
even more accurate determination of the atomic Rabi fre- 
quencies: the outer sideband peaks' distance is given by 
D s (0) = 4G = + g\ x (in the secular limit [25] and 

when the x-ray laser is on resonance, i.e. its detuning 
A = w 3 i — uj x from the transition frequency is 0). In this 
formula, the optical Rabi frequency <? 2 i is usually known, 
therefore, determining its counterpart g 3 \ for the x-ray 
transition is only limited by the accuracy of measuring 
the peak distance D. Fig. [2] (b) shows the ratio of the 
width of the narrowed outer lines to their distance D s (0) . 
As shown, this ratio, characterizing the relative accuracy 
for the determination of atomic multipole moments, can 




821 [eV] 

FIG. 2: (Color online.) (a) Density plot of the fluores- 
cence spectrum (logarithmic scale, arb. units) as a function 
of the fluorescence photon frequency ujf with respect to the 
x-ray transition frequency ^31 (abscissa) and the laser de- 
tuning A = uj x — uizi (ordinate), with the frequencies nor- 
malized by the T31 rate. The parameters are for Bi as in 
Fig. [l] (b) Continuous (blue) curve: ratio of the interference- 
narrowed width Tsb of the outer sideband peaks to their dis- 
tance -Ds(O) = AG — A*J gf-L + g^x as a function of the optical 
Rabi frequency <?2i , with further parameters for the Bi three- 
level system as given in the third line of Tablell] Dashed (red) 
curve: deviation of the sideband distance D S (A = I^i) from 
D a (0). Dotted (green) curve: deviation of the exact sideband 
distance D(0) from its value in the secular limit D s (0). 

be improved by several orders of magnitude for higher 
optical laser intensities. 

Our calculation, in analogy with the derivation of 
Ref. [TT] for zero detuning, shows that the outer side- 
band distance for non-zero detuning is given by D S (A) — 
4G + f (AR - 3i? 2 ) (A/G) 2 + O ((A/G) 4 ). This weak 
dependence is also illustrated in Fig. [2] (a). Hence, the 
experimental sensitivity on the potentially inaccurately 
known detuning may be reduced by orders of magnitude 
by increasing the optical intensity (Rabi frequency), as 
also shown in Fig. [2] (b). The multipole matrix elements 
of the ionic transitions can thus be determined in princi- 
ple to a high accuracy on the order of 10 -4 -10 -6 , once 
the intensity of the driving field is accurately known. 
Conversely, the intensity may be measured to high ac- 
curacy if the multipole moments are reliably known from 
an independent experiment (e.g. lifetime measurements). 
At the same time, knowing the Rabi frequencies, the de- 
pendence of sideband positions on the x-ray detuning 
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TABLE I: Parameters for 2s-f->2p 3 / 2 transitions in the Li- 
like ions 203 T1 78+ (o; 3 i=2236.5 eV), 209 Bi 80+ (2788.1 eV) and 
235tj89+ (4459 4 e v). Optical transition energies (w 2 i, [23]). 
natural line widths (r3i, I^i) and Rabi frequencies (g3i, <? 2 i) 
as well as the interference-narrowed outer sideband width 
(Tsb) °f the x-ray transition are given for the laser inten- 
sities I B , I . x(y) stands for x x 1Q V . 
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could allow to measure in principle the ionic x-ray tran- 
sition energy in an independent way. 

Our above results have been demonstrated on the ex- 
ample of HCI with non-zero nuclear spins, i.e. when hy- 
pcrfine splitting of the ground state occurs. However, 
certainly, the results may be applied to further three-level 
systems. For example, such configurations may also be 
prepared by applying (strong) external magnetic fields, 
which gives rise to a large Zeeman splitting of the ground- 
state level, addressable by long-wavelength coherent ra- 
diation such as masers (or even CO2 lasers). Further- 
more, the results can be generalized to other physical 
systems with high transition energies, such as electro- 
magnetic transitions in nuclei. In this setting, nuclear 
multipole moments and transition energies may in prin- 
ciple be determined by an independent method. 

Laser systems with photon energies of up to a few keV 
(in the range of Li-like transitions) are presently avail- 
able [57] , allowing to excite elements as heavy as U 
and observing emission lines of sub-natural linewidths. 
Upcoming laser facilities are expected to increase the fre- 
quency limit to the order of tens of keVs (e.g. [IQ]), per- 
mitting to directly address the most relativistic heavy 
H-like systems or even nuclear transitions. Table [I] lists 
values for some elements and atomic transitions. 

In summary, a fully relativistic ab initio theory of the 
bound dynamics of atomic systems in laser fields rang- 
ing to the x-ray domain has been developed. The bare 
atomic states are constructed from solutions of the Dirac 
equation. This approach allows for exploiting the sensi- 
tivity of inner-shell electrons to relativistic electron corre- 
lation, QED and nuclear effects in strong Coulomb fields. 
As a demonstrative example, a means to determine ionic 
transition multipole moments and frequencies via a three- 
level configuration driven by an x-ray and an optical field 
has been put forward. Current or near-future laser sys- 
tems are expected to increase the accuracy of multipole 
moment determinations from the current 10 -3 level (via 



lifetime measurements) to the 10 -4 range or better. Fur- 
thermore, the undesirable trapping of atomic population 
in a long-lived metastable state - naturally occurring 
in certain three-level systems - can be reversed by the 
scheme presented here. Other scenarios developed for 
the quantum control of non-relativistic resonance fluo- 
rescence emission [T2lfT§] are anticipated to yield further 
improvement of detection and accuracy. 
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